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Abstract 

We show that models of inflection point inflation exhibit a phase transition from 
a region in parameter space where they are of large field type to a region where they 
are of small field type. The phase transition is between a universal behavior, with 
respect to the initial condition, at the large field region and non-universal behavior 
at the small field region. The order parameter is the number of e-foldings. We find 
integer critical exponents at the transition between the two phases. 
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1 Introduction 



In this paper we present and study a phase transition (PT) phenomenon in models of 
cosmic inflation. The PT takes place when we interpolate in a specific way between 
models of large field inflation and small field inflation. 

We shall use two definitions to make the distinction between large and small field 
models and refer to them later as first and second categorizations: 

I. Potentials in which the slow roll conditions for the canonically normalized inflaton 
field are satisfied over a distance in field space which is large compared to the four- 
dimensional Planck mass scale are called "large field" potentials, while potentials 
with only a small region where the slow roll conditions are met are "small field" 
potentials. 

II. Adopting the categorization first introduced in [I], we distinguish between potentials 
according to their curvature, V", at the slow roll region. Potentials satisfying V" > 
are called "large field" potentials, while "small field" models have V" < and "linear 
field" models satisfy V" = 0. 

These two categorizations, though in some cases complementary, are not equivalent. For 
example, the recent model [2] involving a potential is considered a large field model 
according to the first categorization, but a small field model according to the second. 
Similarly the even more recent model of [3] is a linear field model according to the second 
categorization and a large field model according to the first. Nonetheless we shall make 
use of both categorizations later on. 

Models of large and small field inflation are quite different: 

• For single field models of inflation, the tensor-to-scalar ratio r and the spectral index 
n s are related by: 

8 „ \ 16 V" 
r =3 (1 "^ ) + W (L1) 
This relation shows that positive curvature models, i.e. large field models with 
V" > are expected to generate larger values for r and can be constrained by 
its measurement through gravity wave experiments, while negative curvature small 
field models predict small values for r. This implies that experimentally models 
with V" > are more severely constrained pi]. 
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• In large field models, where the slow roll region is extended, the initial conditions of 
the inflaton generically do not play an important roll (the field will either start in 
slow roll conditions or slow down enough in the slow roll region to allow for enough 
e-foldings to accumulate and since inflation erases all memory of what happened 
before the last e-foldings, the exact initial conditions play no important role). In 
small field models, however, unless the inflaton starts out near the slow roll region 
of the potential, it will overshoot it without inflating the universe. Simply put, the 
Hubble friction does not have enough time to slow down the inflaton at the slow 
roll region. 

• Constructing large field models in string theory is quite difficult (see, however, 
[21 [3]). Small field models are generally easier to construct, largely since the slow 
roll conditions can be satisfied accidently in a small region. 

Since the two types of models are distinctly different, one might ask if the transition 
between them is smooth or not. Namely, is there a PT between the two? To explore this 
possibility we have to study the physics along a trajectory which interpolates between 
small field and large field inflation. The problem is that the space of all possible models of 
single field inflation has an infinite dimension, and so there are infinitely many trajectories 
from the region of small field inflation to the region of large field inflation. Thus to proceed 
in a meaningful way we have to consider a slice of this space that truncates the number 
of dimensions. In the present paper we consider the simplest possibility of doing so via a 
one dimensional slice which is defined in the following way. We rescale <fi in the potential 
V{<f>) 

0^A0, (1.2) 

while keeping the kinetic term canonically normalized. From the definition of the slow-roll 
parameters 

1 fV'V V" 

^Av) • ' = t- (L3) 

and our first categorization above, we see that for any non-singular potential, A — > oo 
gives a small field model while A — > yields a large field model. 

The total number of e-foldings obtained at the slow roll region (the region where e < 1) 
is denoted by N. As we vary A we expect to find the two extreme values 

lim N{\) = 0, limiV(A) = oo. (1.4) 

A^oo A^O 
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This suggests that it is natural to think of N(X) as an order parameter between the two 
phases of large and small field inflation. 

The main point of the paper is that contrary to what might be expected, in some 
cases, in which V(4>) is regular and has an inflection point, there is a PT at some finite 
value of A where N suddenly becomes infinite. This discontinuous transition happens 
when the parameter is varied continuously, reminiscent of phase transitions of the second 
kind in Landau- Ginzburg theory. Moreover, we find a scaling behavior near the critical 
point of the form: 

N W ~ (a -\ c y ' A < Ac • (L5) 

As is often the case with scaling behavior, ( 11. 5ft is universal. The sense in which it is 
universal here is that it does not depend on the initial condition of the inflaton. Namely, 
near A c any non-singular initial condition yields (11. 5p . This is due to the attractor property 
of inflation^ This scaling obeys an integer exponent scaling law as we later show and 
explain in detail. 

The paper is organized as follows. In section 2 we demonstrate a PT in models of 
inflection point inflation (IPI). Such models have recently become popular in the context 
of inflation in string theory El EH EU H21 1131 M> EE] and the MSSM [161 El EE] . We 
show numerically that IPI models satisfy the scaling behavior described in ( 11.51) with 
7 = 1. In Section 3 we discuss generalizations of the PT phenomenon. We show that it 
extends to models with an approximate inflection point and check the influence of higher 
powers of (f) n in the potential. For example, it turns out that the critical exponent 7 of 
the scaling behavior directly depends on the lowest degree of non-vanishing derivative in 
the extremal point of the potential around which inflation occurs. Section 4 includes an 
analytic treatment of the PT properties described in Sections 2 and 3. Finally, in Section 
5 we illustrate that a somewhat different scaling behavior also takes place in a scenario 
with a time-independent IPI potential supplemented by a time-dependent contribution 
which helps solve the initial condition problem [12] , This model has a similar PT (also 
with integer scaling behavior near the critical PT point) that occurs at a critical value of 
the initial condition of the time-dependent contribution to the potential. We conclude in 
section 6. 

1 For a general discussion on the attractor mechanism see [5] . The efficiency of the attractor mechanism 
in inflection point inflation models was discussed recently in [6l [7] . 
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2 Phase Transition in Inflection Point Inflation 



In this section we consider the basic example of a PT between models of small and large 
field inflation. The example is of inflection point inflation (IPI) where the inflaton's 
potential can be approximated near the inflection point, at = 0, by 

V{ct>) = V {l + f3^). (2.1) 

The fact that these kind of models are special follows from the second categorization above. 
According to this categorization, for any (3 the models are of large field for sign(/5)0 > 
and small field for sign(/5)0 < 0. Indeed, as we shall see, what determines the transition 
is whether the inflaton spends most of the time at positive or negative values. 

From the point of view of our first categorization above, what determines if the po- 
tential is of large or small field form is the value of (3, which controls the shallowness of 
the potential near the inflection point. For (3 « 1, the potential is of large field form 
(it is very flat near the inflection point and so the slow roll conditions are satisfied over a 
large range in <p space) while for f3 » 1 it has small field form. 

We wish to demonstrate the PT property numerically, by finding a finite A c for which 
N(X C ) = oo. From (12.11) . we see that taking — > A0 and keeping the kinetic term fixed is 
equivalent to keeping <p fixed and transforming 

(3 — ► A 3 /3. (2.2) 

Therefore, to demonstrate a PT we need to find a finite (3 C . 
The equations of motion are 

3H 2 = 1 -^ + V (l + (3<b 3 ), 

4> = -3H<j)-3V f3(f) 2 , (2.3) 

where as usual H = a/a. The classical equations of motion are invariant undei0 

V -> a 2 V , t -> t/a. (2.4) 

This transformation does not change the region (in field space) of slow roll and it leaves 
our order parameter, N — J Hdt, invariant. Thus it can be used to fix Vq = 1. 

2 This symmetry is broken at the quantum mechanical level by the Planck scale. 
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Table 1: The amount of e-foldings generated for different initial conditions (assuming 
(fiinit = 0) for different values of the free parameter (3 (holding Vq = 1). We see that below 
the critical value f3 c ~ 0.7744 inflation is infinite for all initial conditions. 

We are left with the task of solving (12.31) . with Vq = 1, numerically and calculating 
N (while the inflaton is in the slow-roll region) as a function of (3 for different values of 
the initial condition. These results are summarized in Table 1 which shows that there is 
a critical point around (3 C ~ 0.7744 below which N goes to infinity independently of the 
initial condition. 

To make contact with critical behavior and second order phase transitions it is more 
convenient to work with ~ (rather than N) as the order parameter. The reason is that it 
changes continuously from zero to finite values, as shown in figure 1. Indeed this figure 
looks very familiar from discussions on second order phase transitions. We can go even 
further in improving the definition of the order parameter. As we approach the critical 
point, most of the e-foldings are obtained near the inflection point. Since in this limit 
H is roughly constant near the inflection point what determines the behavior of i is the 
amount of time the inflaton spends near the inflection point. This in turn is determined 
by the inflaton's velocity at the inflection point 

I ~ 0(0 = 0) = O . (2.5) 

Therefore, it makes sense to work with 0o as the order parameter. Moreover, since 0o is 
a local operator (while 1/N is a non-local operator) it is easier to numerically calculate 
0o . In the remainder of the paper we shall use these alternatives as the order parameters 
for our analysis. 

In Figure 2 we present the scaling behavior near (3 C both for 1/N and for 0o- We see 
that both cases exhibit scaling behavior when j3 c is approached with excellent accuracy. In 
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Figure 1: Alternative for the order parameter of our PT. The behavior of ■h as a function 
of (3 is shown. The transition at the critical point is continuous, dropping to zero when it 
is crossed and inflation becomes infinite. 

addition, we see that the slope of all graphs shown is extremely close to 1, which implies 
the relation 

0o -^oc^-/^) 1 . (2.6) 

The fact that the critical exponent is an integer strongly suggests an analytic explanation. 
Such an explanation will be discussed in section 4. Before we turn to this explanation we 
consider in the next section some generalizations of this basic example. 

3 Irrelevant and Relevant Deformations 

In this section we consider some relevant and irrelevant deformations of the potential 

flZH). 

We start with the irrelevant deformations. Consider the following potential. 

V{<j)) = V ( 1 + /?0 3 + a ^ ni ) ■ 

V rii>3 ) 
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Figure 2: Scaling of i and 0o with distance of /? from its critical value. Both cases obey 
a scaling law with critical exponent 1. 

As before, we would like to study the PT as we rescale <j) according to eq. (jl.2p . This 
rescaling is equivalent to taking 

(3 -»■ A 3 /3, «i -»■ \ n *a t . (3.2) 

Our numerical simulations show that, as expected, the higher order terms are indeed 
irrelevant in the sense that they do not influence the existence of the critical point and 
the value of the critical exponent. 

Only when the "marginal term" vanishes (J3 = 0) the critical exponent is not 1. When 
the leading term in the potential is 

(l + «0 2n+1 ), n = l,2,... (3.3) 

the critical exponent is In — 1. Potentials in which the lowest non-zero derivative is even, 
i.e. of order 2n, exhibit no critical behavior at all. We shall return to this in the analytic 
discussion in section 4. 

Next we consider a relevant deformation of the form 

V = V (l + 6(f) + (34> 3 ). (3.4) 
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Figure 3: Scaling of ■h with distance of parameter from its critical value. Left: Ap- 
proximated IPI potential with non-zero fifth derivative, critical exponent is 3. Right: 
Approximated IPI potential with non-zero seventh derivative, critical exponent is 5. 



This time the rescaling (jl.2B is equivalent to 

(3 -> X 3 (3, 5 -> X6. (3.5) 

It is easy to see that for any 5^0 there is no PT of the kind discussed before. Namely, 
there is no finite A c for which iV(A c ) = oo, and hence no scaling behavior of the form 
(ll.5j) . From the point of view of 0o as the order parameter this is clear since the minimal 
value of 0o is obtained in the slow roll approximation and is linear in 5. Therefore, we 
are never in the phase where the order parameter vanishes. 

This can also be seen when considering 1/N as the order parameter. The maximum 
number of e-foldings, obtained by tuning the initial condition, is [9] 



Nm ~ = n yw (3 - 6) 

Again we see that for any finite value of S the order parameter, 1/N, never vanishes. 

Still, as is often the case in Landau- Ginzburg theory, some aspects of the PT remain. 
Here we find that for any 5^0 there is a A c above which N depends on the initial 
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condition, and to maximize iV we have to fine tune the initial condition. On the other 
hand for A < A c any non-singular initial condition gives N = N max . From ( 13. 5ft and ( 13.61) 
we find the universal (with respect to the initial condition) behavior for A < A c : 

N{\) = N(X C ) 02 V . (3.7) 

The fact that N(0) = oo is trivial. What is interesting here is the fact that ( 13.71) does not 
depend on the initial condition, and that it scales like 1/A 2 . 



4 Analytic Discussion 

In this section, our aim is to find an analytic explanation for the integer critical exponents 
we have found numerically. In order to do that, we have to simplify the equation of motion 
keeping only the most relevant features. Since near the critical point the inflaton spends 
most of the time near the inflection point, This means that as j3 approaches (3 C we can 
approximate H by a constant. This assumption allows us to replace our order parameter 
with - where t is the time spent by the inflaton in the vicinity of the inflection point 

(since N e _f Q i ds = f Hdt = H(t end — t start )). We have also seen that irrelevant terms ( 

t start 

<fr n terms with n > 3) are indeed irrelevant, in the sense that adding them to the potential 
does not modify the critical exponents. Therefore we expect to find a similar PT with 
the same critical exponents for the following equation 

4> = -3H <f) - 3VW 2 , with 3# 2 = V . (4.1) 

This equation approximates (I2.3P near the inflection point and is invariant under the same 
transformation (12.41) . Again we use this symmetry to fix Vq (or H ). This time we find it 
useful to fix Ho = 1 /3 which leaves us with 

0=-0-/30 2 . (4-2) 

Unlike the original equation, fl 2 . 3 j) . this equation is also invariant under 

<j> -> C(f), (3^ (3/C 2 . (4.3) 

This means that solutions to (14.21) with the initial condition <p(t = 0) = 0o an d 0(t = 
0) = are parameterized by the combination <J)q(3 which is invariant under (14.31) . Thus 
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to describe the PT we can either fix </>o and vary (3 or fix (3 and vary 0o- We find it most 
convenient to fix (3 = 1 and vary the initial condition 0o- 

Our goal in the rest of this section is therefore to show that the equation 



> 2 . 



(4.4) 



exhibits a PT as a function of 0o and to find the relevant critical exponent. To meet this 
goal we use the phase-space approachjfl We rename the axes as x = 0, y = and reduces 
the second order differential equation to a system of two first order ones 



x = y, 

2 



y = -y-x . (4.5) 



By examining the derivative 



dy = jf = -y-x 2 

dx £j| y 



(4.6) 



we deduce a few general conclusions: 

• When crossing the x axis (y = 0) the slope is always infinite. 

• When crossing the y axis (x = 0) the slope is always —1. 

• When crossing the contour y = —x 2 , the slope is always zero (except at x — 0). 

• For negative y, outside the contour y = —x 2 , the derivative is positive and so the 
trajectories are monotonically increasing. 

• For negative y, inside the contour y = —x 2 , the derivative is negative and so the 
trajectories are monotonically decreasing. 

This explains the trajectories in Figure 4, from which we see that indeed there is a critical 
initial condition point and that orbits starting at different points on the x-axis (i.e. initial 
conditions with zero velocity) either converge to the stable point at x = y = or overshoot 
it and asymptotically approach x,y — > — oo . 

These general considerations, however, are not sufficient to determine the critical 
exponent. To derive the critical exponent we have to analyze the behavior near x — 0. 



3 We are grateful to Barak Kol for suggesting this to us and for deriving some of the results below. 
This approach was used recently in [7j. 
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Phase Space Digram for: x+jrt-x^O 
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Figure 4: Phase Space trajectories of solutions to (}4.4j) for different initial conditions. 



The region near x = can be approached in two different ways: 

(I) We can neglect the x 2 term. Then the second order equation becomes linear, with 
the solutions 

x = cie~* + c 2 . (4.7) 

This corresponds in the original problem to keeping the Hubble friction term while ne- 
glecting the slope of the potential. Namely, we are in the de-Sitter approximation. This 
solution is valid to the right of x = and we see that for c 2 = the trajectory will reach 
the origin x = y = 0ait — > oo, coninciding with the critical trajectory beginning exactly 
at the critical point on the x axis. 

(II) We can also neglect the x term. Then the equation becomes first order and we 
have the solution 

x = -±-. (4.8) 
t + c 

This corresponds in the original problem to taking the slow roll limit. This solution with 
c = describes a trajectory that starts at the origin at t — > — oo and moves to the left 
of x = 0, coinciding with the curve along which the derivative vanishes (the curve 
y = —x 2 with x < 0). 

Now, in order to find the critical exponent, we need to consider the behavior of a 
trajectory that starts at some small distance e' from the critical point on the x axis. Such 
a trajectory will never reach the origin, but will eventually cross the y axis below zero. 
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However, when it approaches x — 0, the critical trajectory is a good approximation for it 
(see Figure 5). Therefore, we start with solution (I) 

x = ce~ l - e , y — -ce - *, (4.9) 

parametrizing the distance from the c2 = solution by e. Here it is important to note 
that since the equations are regular in this region, the relation between e' and e is linear. 
Therefore, finding a relation between t near the x = and e should give us the critical 
exponent. The trajectory will cross the x = point at y — — e. We shall continue with this 
solution until it becomes unsatisfactory, at which point we will replace it with the second 
solution. The approximation (I) clearly breaks down when x ~ x 2 ~ e 2 , which happens 
at 

t* ~ log (4) , x* ~ e 2 - e ~ -e , ~ -e 2 . (4.10) 

Thus the first approximation is reliable up to x ~ — e, which, as is clear from Figure 5, 
is indeed in the vicinity of the curve along which approximation (II) is valid. However, 
as we saw above, the solution trajectory crosses this curve parallel to the y axis and 
therefore must deviate from it (mush as it deviated from the critical trajectory to the 
right of x — 0). We can therefore continue with the second approximation (I), starting at 
x ~ — e, until it, too, becomes invalid. To find when this occurs, we start with 

x ~ —e — 5x , x ~ — e 2 , (4-H) 

and stop when x ~ x ~ e 2 , at which stage neglecting the x term is no longer a good 
approximation. We find 

„ = dydx = -y-x 2 y = _ y _ x 2 =4 2_ f _ 26xe _ {6x) 2 „ 6xe (412) 
ax at U 

which is no longer negligible once 5x ~ e or x ~ — 2e. 

We conclude that the time spent by the inflaton in the vicinity of the inflection point can 
be divided into regions with different characteristics: in the interval to the right of zero 
it behaves as t rig ht ~ — In e, while in the interval [— 2e, — e] it satisfies 

Ueft ~ -■ (4.13) 

Since for e << 1 we have | lne| << ^, this implies that indeed the critical exponent in 
the relation between | and (xmu — x cr u) is 7 = 1. 
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Figure 5: Phase Space trajectories of solutions to ( 14.41) for different initial conditions. The 
curve y = —x 2 , along which the derivative ^ = is also plotted. 

We end the section by considering very briefly the higher order cases of section 3. 
Using the same considerations as before we find that a potential of the form (13.31) leads 
to 

x + x + x 2n = 0. (4.14) 

Drawing it in phase space (see Figure 6) we see the same qualitative behavior as for the 
previous potential. Repeating the approximations we did in the last section (solution 
(I) will obviously remain the same, while in solution (II) we will have x ~ fin-i ) ^ ne 
logarithmic behavior to the right will remain, but to the left we will get the relation 
x ~ t 2n-i which indeed explains a scaling law with an integer exponent of 2n — 1. 
On the other hand, for potentials with <fr 2n as the lowest non-zero terms we get 

x + x + x 2 " 1 ' 1 =0. (4.15) 

Note that this equation, unlike (14.141) . is invariant under x — > —x. This implies that there 
is no PT. In Figure 7 this is illustrated for n = 3 where we can see that no initial condition 
leads to a critical solution. 
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Phase Space Diagram for: x 4- x 4- x =0 



Phase Space Diagram for: x + x + x =0 




Phase Space trajectories 
1 Trajectories for some specific initial conditions 
Contour along which dy/dx=0 



Phase Space trajectories 
■ Trajectories for some specific initial conditions 
Contour along which dy/dx=0 



Figure 6: Phase Space trajectories of solutions to (I4.15P and (14.141) for different initial 
conditions. 



5 PT in a Time-Dependent IPI Model 

So far we have discussed the appearance of phase transitions in models of time- independent 
IPI. In this section we show that this phenomenon extends to cases of time-dependent 
potentials, although in this case the critical parameter is not a parameter in the potential 
or the initial condition for the inflaton, as we shall momentarily see. To demonstrate this, 
we focus on a stringy inspired time-dependent setup that was introduced in [12]. The 
relevant equations of motion are 



3H 2 



n 



\4> 2 + v static + \ ; 

-3Hn, 



dynamic} 



(5.1) 



+ 3H<j) 



static 



+ v c 



dynamic) j 



where V sta tic is a small field potential (according to the first categorization) with an in- 
flection point. 

As is clear from the discussion above, the problem with such a V s tatic is that for a 
generic initial condition the inflaton will overshoot the inflection point and the universe 
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will not innate. The resolution proposed in (12] to the proble nfl is that particles with 
mass that depends on the expectation value of the inflaton, and satisfies 

dVstatic dM 

induce a time-dependent potential 

Vdynamic = Tl(t)M(<p) (5.3) 

that slows down the inflaton at the inflection point and resolves the overshoot problem 
(see Figure 7 for more details). 

In [12] it was observed that the initial value of Uq needed in order to generate a large 
number of e-foldings is quite small. Here we make the point that there is a critical value 
of Uq which we denote by n c , above which N = oo. In addition, we find that as n 
approaches n c from below we obtain a scaling behavior 

N , (5.4) 

n c - n 

which implies that in this case the critical exponent 1 is also obtained with respect to the 
initial condition. 

Before discussing the numerical results, let us present some of the technical details. 
The potentials we worked with are the ones of [T2] that are typical in stringy compactifi- 
cations. 

V static = a>i exp(-jia0) + a 2 exp(-j 2 a0) + a 3 exp(-j 3 a</>) 

1 



ji = 12, j2 = 10, js = 8, a 



24 

^inflection? r ^inflection ) ) 

(5.5) 



Oi 2 4 02 _ u 2 

q inflection' r inflection! 

Q>3 "J Q>3 3 



where ^inflection — e , and L is the size of the compact manifold. The particles with 

masses that depend on the inflaton and satisfy (15.21) are branes that wrap the compact 
manifold. They induce 

Vdynamic = n(t) exp(3«0). (5.6) 

(We took the initial condition to be <fi = <fi = 0, which corresponds to L = 1, L — 0.) 



Similar resolutions were proposed in the context of the overshoot problem associated with moduli 
stabilization [19] in l20l. 
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Figure 7: A heuristic demonstration of how a time dependent potential, that scales 
like l/a 3 (t), can resolve the overshoot problem of IPI. The static piece of the potential 
is denoted by the red/solid line and the time dependent piece of the potential by the 
green/dashed line. In (a),(b) the static potential is steep, but since there has not been 
significant expansion the time-dependant potential is able to balance it, preventing <p(t) 
from acquiring a large velocity. Heuristically, 4>(t) (denoted by the blue/filled circle) fol- 
lows the dilution of the time dependent potential. In (c) <fi(t) enters the shallow region 
(where the slow roll parameters are small) where inflation takes place and the time de- 
pendent potential starts to slow down exponentially. In (d) 4>(t) is dominated solely by 
the static potential but its velocity is now low enough to allow prolonged inflation. 
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Table 2: N as a function of n and L in fi ection . We see that for any L in fi ection there is a 
critical value of no above which N = oo. 

Now, solving the equations above numerically, we find a critical initial condition for 
n(t), above which the inflaton is stuck at the inflection point (see table below). Near the 
critical point, observing t (the time spent near the inflection point) instead of N, as we 
did in Section 4, we notice the same scaling behavior, obeying a critical exponent of 1. 
This is plotted in Figure 8. Our simulation shows that the critical exponent does not 
depend on the details of V sta tic and Vdynamic- For example, replacing the 3 in (15.61) by, say, 
2 would change n c but will not change the critical exponent. 
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Scaling of t as a function of log(n-n ) 

7.4 
7.3 
7.2 
7.1 
4J 7 
6.9 
6.8 
6.7 
6.6 

-14. „ . 

log(n-n ) 

c 

Figure 8: Scaling behavior of the string Theory model near the critical point in the initial 
condition of the time-dependent contribution to the potential. Plotting the time spent in 
vicinity of inflection point vs. the Log of distance from critical point, we see the same 
scaling behavior encountered before, obeying a critical exponent close to 1. 

6 Conclusions 

The goal of the paper was to show that in IPI models, that have become quite popular in 
string theory recently, the efficiency of the attractor mechanism varies in a discontinuous 
fashion as we vary various parameters continuously (as happens in phase transitions of 
the second king in Landau-Ginzburg theory). The parameter space of these models can 
be divided into two regions. In one (the small field region), N depends on the initial 
condition and generally speaking, to obtain a large N the initial condition should be fine 
tuned. In the large field region the efficiency of the attractor mechanism is maximized 
and all non-singular initial conditions yield the same N. The surprising aspects of our 
findings are: 

(I) The transition between the two regions happens at finite values of the parameters that 
define the potential. 

(II) In some cases critical behavior takes place at the transition between the two regions, 
obeying an integer exponent scaling law, which is shown and explained above. 

These findings shed new light on the approach one should take in examining an IPI 
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model in terms of the required fine tuning of the parameters of the model and the initial 
conditions. A good example is the time-dependent string theory model, where it is now 
evident that less care in fine tuning of the initial density of the stringy particles is needed 
in order for the proposed mechanism to work, since below a specific (and quite low) critical 
value, all values are equally sufficient. 
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